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UNIFORMIZATION OF NONLINEAR HAMILTONIAN 
SYSTEMS OF VLASOV AND HARTREE TYPE 

V.P. BELAVKIN AND V.P. MASLOV 


Abstract. Nonlinear Hamiltonian systems describing the abstract Vlasov 
and Hartree equations are considered in the framework of algebraic Pois- 
sonian theory. The concept of uniformization is introduced; it gener¬ 
alizes the method of second quantization of classical systems to arbi¬ 
trary Hamiltonian (Lie-Jordan, or Poissonian) algebras, in particular to 
the algebras of operators in indefinite spaces. A functional calculus is 
developed for the uniformized observables, which unifies the calculi of 
generating functionals for classical, Bosonian and Fermionian multipar¬ 
ticle variables. The nonlinear kinetic equation of Vlasov and Hartree 
type are derived for both classical and quantum multiparticle Hamil¬ 
tonian systems in the mean field approximation. The possibilities for 
finding approximate solutions of these equations from the solutions of 
uniformized linear equations are investigated. 


1. Introduction 

Recently, much interest T_ has been devoted to studying the region of 
the spectrum of the Hamiltonians of certain self-interacting boson fields 
corresponding to particle-like (soliton) solutions of the classical mean field 
non-linear equations. Such excitations can be naturally interpreted as com¬ 
posite stable quantum particles that have finite energy even in the classical 
limit. As a comparison |2J of the results for exactly solvable models showed, 
there is a direct correspondence between the quantum [3. and classical [3] 
particle-like solutions. Frequently, it is simplest to obtain such solutions di¬ 
rectly for the quantum case, and not for the corresponding classical case, for 
which they were obtained only very recently by the solution of the inverse 
scattering problem. 

This correspondence can be viewed as a heuristic basis of an asymptotic 
method for finding soliton solutions of nonlinear equations in not only classi¬ 
cal held theory but also for the difference equations that describe the dynam¬ 
ics of nonlinear crystals, and also certain nonlinear equations of quantum 
electrodynamics. This method, which we call the uniformization method, 
can be regarded mathematically as a new method of linearization leading to 
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asymptotically exact solutions of nonlinear equations in the large number 
limit by the infinite increase in the number of degrees of freedom. 

Applied to ordinary nonlinear Schrodinger equations, the uniformization 
method essentially reduces to second quantization and the construction of 
a quasiclassical soliton solution from the solution of the n-particle equation 
for sufficiently large n. However, this interpretation no longer holds for the 
spaces with indefinite metric. Applying this method to the nonlinear equa¬ 
tion of quantum electrodynamics ([Sj, p.ll) in which the local time enters 
as a fourth coordinate, we arrive at a many-time formalism that reduces 
the solution of these nonlinear equations to the solution of ultrahyperbolic 
equations that may not have independent physical meaning. 

Therefore, we set forth the uniformization method on an abstract alge¬ 
braic basis without referring to its physical interpretation which is contained 
in the concrete functional realizations. As the object of uniformization, we 
choose the abstract Vlasov equation, which generalizes the classical and 
quantized Vlasov equations considered in 0, and also equations that re¬ 
duce to equations of Hartree type in spaces with indefinite metric. 

For the classical case, the uniformization method is essentially the tran¬ 
sition from a field system described by the Vlasov kinetic equation to a 
system with variable number of particles, i.e., it is the procedure which is 
the inverse of the derivation of the Vlasov equation from the Bogolyubov 
chain of n-particle equations In the quantum case, the uniformization 
method gives a unified algebraic basis for the method of second quantization 
of bosons and fermions with a superselection rule that avoids the consider¬ 
ation of unphysical - odd - observables that do not commute with the par¬ 
ticle number operator. The developed calculus of representing functionals 
for such even operators combines and replaces the calculus of symbols of 
commuting variables (Bargmann j^) and anticommuting variables (Berezin 
jjj). In such an approach, the difference between bosons and fermions does 
not reside in the operators but only in the multiparticle tensor-states. 

In the general case it is shown that the algebra of polynomial observables 
of the uniformized system is isomorphic to the direct sum of the tensor 
powers of the Hamiltonian algebra generated by the elementary observables 
of the uniformized system. In the case of algebras of operators in a space with 
indefinite metric, such a representation generalizes the Fock to a pseudo-Fock 
representation. 

First, we set forth the basic concepts and the definitions relating to non¬ 
linear Hamiltonian systems of Vlasov and Hartree type, and also the uni¬ 
formization method, showing that such nonlinear kinetic equations naturally 
appear as characteristic equations in the mean-field approximation of the 
uniformized Hamiltonian infinite systems. Then we consider the possibility 
of finding approximate solutions of these nonlinear equations, in particu¬ 
lar approximate soliton solutions from the corresponding solutions of the 
uniformized equations. 
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2. Abstract Vlasov and Hartree Equations 

2.1. We consider a dynamical system whose states are described by the 
elements p of a real linear space C as the predual space of linear functionals 
A i—> (p, A) on an abstract Hamiltonian algebra A with respect to some 
bilinear form (p, A) having real values for all p £ C and A = A* 6 A 1 . By 
the Hamiltonian algebra we mean any complex-linear unital *-algebra A in 
which besides the “ordinary”, not necessarily commutative and associative, 
product A ■ A* £ A, there is also defined the Poisson brackets {H, A} as a 
bilinear antisymmetric map A x A —> A with the property of continuous 
derivation 

{H, A • A*} = A ■ {H, A*} + {H, A} • A* 

having anti-Hermitian values {A*, A }* = — {A*, A} in contrast to {A ■ A*)* = 
A • A* (for more details see 7|). If the ordinary product is associative, the 
Hamiltonian algebra is said to be Poissonian. In the commutative case 
A - A* = A* ■ A the real part 3?A = {A = A*} of the Hamiltonian algebra is 
a real Hamiltonian subalgebra such that A can be considered as a complex- 
ification of 1Z = KA which is usually taken to be a real Lie-Jordan algebra 
fSj. If the algebra A is both commutative and associative with respect to 
the product A • B like the algebra of smooth scalar-valued functions on a 
Hamiltonian phase space, it is called classical Hamiltonian, or, equivalently 
Poissonian algebra. 

We give the name Vlasov Hamiltonian to every real element H = H * £ A 
that may depend continuously on p £ C as the derivative 6 j(p) = H(p) 
of some, in general nonlinear A-smooth functional 7 : p 1—> 7 (p), called 
Hamiltonian functional on £ (7 has real number values 7 (p) £ R). By 
A-smooth we mean such smooth functional which has the derivative as a 
continuous functional p 1 —> <Py (p) defined for each p £ C as an element H (p) 
of A such that 

(2.1) (a, 5) 7 (p) := d ^ p g + t ta) | t=0 = (a, H(p )), Va £ C. 

Let us assume that the linear operator p 1—> {/>, H} is well-defined for each 
value H = H (p) of a Vlasov Hamiltonian H as the predual to the derivation 
A} on the algebra A such that 2 

(2.2) ({ p,H},A) = (p,{H,A }), Vp £ £, A £ A. 

We give the name “abstract Vlasov equation” to the evolution equation in 
the state space C 3 p(t) of the form 

(2-3) ^ = 

^One may assume that £ is a Frechet space with the complexified dual of all linear 
continuous functionals p 1 —> (A, p) given by a *-algebra A. 

2 To this end one may assume that the derivation A 1 —> {H, A}, as well as the action 
A >—> H ■ A are a (A, £)-continuous for each H £ A. 
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where H(t,p) is the derivative 8^{t,p) of a nonlinear time-dependent A- 
smooth functional 7 (t, p) which may, in general, depend on time t. 

Proposition 1. Equation \2.A) is a classical Hamiltonian system with re¬ 
spect to the Hamiltonian functional y(f, p) which has at least one linear 
integral of motion u(p) = (p, I) generated by the identity I E A. 

Indeed, for any differentiable functional a : p 1 —> a(p) that does not 
depend explicitly on t we have in accordance with (fOll 

(2.4) = ({p(t),H(t,p(t))} ,Sa(p(i))) . 

Using the definitions Q of the dual operator A 1 —> {H, A} and remember¬ 
ing that H = 8 r y(p), we can write (12.41) in the form ^ = { 7 , a} d , where 

(2.5) { 7 , a} d (p) = (p, {8r/(p),6a(p)}) 

is a classical Poisson bracket with respect to the pointwise product a/3 since 

{ 7 , (p) = a (p) { 7 , (3} d (p) + { 7 , a} cl (/ o)(3 (p). 

This means that the algebra 21 c i of smooth functionals a : C —► C is a classi¬ 
cal Hamiltonian algebra with commutative and associative Jordan product 
(a • (3) (p) = a{p)(3{p) and Lie product (12.51) which generates on the man¬ 
ifold C a symplectic structure with respect to which ifO) is simply the 
Hamiltonian equation with the Hamilton functional 7 (t,p). Substituting 
v(p) = (p,I) as a (p ) i n (12.51) and remembering that in accordance with the 
derivation property we have {H, /} = 0 for any H £ A, we find that u(p) is 
an integral for (12.51) : { 7 , u} d = 0 . 

Note that the elements a of the algebra 2t c / can be interpreted as nice 
(smooth) observables of the described classical system, taking real or com¬ 
plex values a(p) in the states p E £, which should be regarded as deter¬ 
ministic point-states corresponding to zero variances |a | 2 (p) — \a(p)\ 2 = 0 
simultaneously for all a E 2t c ;. The Hamiltonian algebra A, whose elements 
generate “elementary classical observables” - the linear functional observ¬ 
ables a a = {p, A) - is here not assumed to be classical. 

2.2. Actually to the Vlasov equation itself there corresponds only the case 
when A is the classical Poissonian algebra of smooth complex functions 
A (q, p) on the single-particle phase space M. 2d 3 (q, p) with ordinary point- 
wise multiplication and involution given by complex conjugation A * (q, p) = 
A (q, p), and the classical Poisson brackets { A *, A} = i [A*, A\ d , where 
[A*, A\ d (q, p) = i (<9 q A* • d p A - d Q A ■ d p A *) (q, p) 

is the classical commutator, i.e. the Poisson bracket divided by i = \/—l. 
The space C of real density functions p (q, p) absolutely integrable in the 
product with any A E A is dual to A with respect to the phase space integral 

{ P ' A )= . 1.1 


p (q, p) A (q, p) dqdp. 
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Thus the classical Vlasov equation m can be regarded as the Liouville 
equation with Hamilton function Lf(p,q, t,p) that in general continuously 
depends on the density function p (actually in the original Vlasov equation 
H (p, q, t, p) is an affine function of p corresponding to a quadratic form of 
the functional p i—> 7 (t, p)). 

One can obtain a different class of the classical Hamiltonian systems of 
the form (l2Vl) by considering the quantized Vlasov equations as in 0 , or 
even a more general equations when A is an operator algebra in a complex- 
linear Hilbert or pre-Hilbert space TL with scalar product (-|-), taking C as 
the space of density operators p dual to A with respect to the trace-pairing 
(p,A) = Tr (pA). The Poisson bracket with respect to the noncommutative 
operator product in A can be defined by the usual commutator [A *, A] = 
A*A — AA* as {A *, A} = i [A*, A], so that (12.111) takes the form of a quantized 
Liouville, or von Neumann equation 


dp(t) 

dt 


= i [p(t),H(t,p(t))j , 


with p-dependent Hamilton operator H(t,p). Note that the real part 5L4 of 
the nocommutative operator algebra A is not invariant under the associative 
operator product AB, but it is a nonassociative real Hamiltonian algebra 
with respect to the symmetrized operator product A ■ B = \{AB + BA). 

One can also have intermediate semiquantized cases obtained, for exam¬ 
ple, by partial quantization of the Vlasov equation. In all these cases, if the 
states p satisfy the positivity condition (p, A*A) ^ 0 and the normalization 
condition (p, I) = 1 they can have a probability interpretation by regard¬ 
ing the forms (p, A) not as exact values of elementary observables but as 
expected values of the results of individual measurements of these observ¬ 
ables. However this probabilistic interpretation cannot be extended onto 
the nonlinear kinetic equation even though they preserve the positiv¬ 

ity of p and its normalization, since they do not preserve statistical mixtures 
except in the linear case. As we shall see, the the nonlinear kinetic equa¬ 
tions describe not the individual but a collective dynamics of the infinite 
ensemble of the corresponding classical or quantum particles in the mean 
field approximation. 

Finally, our treatment also encompasses the (general) case when A is an 
operator A-algebra in a space 7i with indefinite metric, i.e., when (•)•) is an 
arbitrary nondegenerate Hermitian form (indefinite scalar product). In this 
case, the set of operators of the form AAA does not form a cone, and the 
states p cannot have a probability interpretation of the individual particles. 
They can only be regarded as pure states of the corresponding classical field 
system, or as a quasi-classical limit of an infinite quantum field system. 
Numerous examples of such a situation arise in a quasiclassical treatment of 
boson and fermion fields with indefinite commutation relations in j 2 j. 
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Suppose ip i—> ip* is a canonical mapping of the space TL into the dual 
space TL * defined by the equation ip* p = (p\ip) , £ 7i, and '0'*/’* are one¬ 

dimensional operators in TL acting in accordance with the formula ipip*p = 
{p\ip) ip. The space C generated by the operators ipip* is the dual to the alge¬ 
bra of self-adjoint operators in TL with respect to the bilinear form ( p , A) gen¬ 
erated by the form (Aip 1^) on the elements p = ipip*. We consider (12.31) with 
the commutator Poisson bracket {■,•} = i and the self-adjoint Vlasov 
Hamiltonian H(t,p) = H(t,p)*. On the manifold M = {ipip* : ip £ TL} C C 
of the one-dimensional operators p = ipip* it takes the form 

(2.6) = H(t,ip(t)ip(t)*)ip{t)ip{t)* -ip(t)ip(t)*H(t,ip(t)ip(t)*), 

which decomposes into a pair of mutually adjoint equations in TL - the 
abstract Hartree equations 

( 2 -7) = H(t,ip(t)ip(t)*)ip(t), -i = ip(t)*H(t,ip{t)ip(t)*). 

This generalizes the Schrodinger equation for the Hartree Hamilton operators 
H(t,ipip*), the Vlasov Hamiltonians on M that depend on one-dimensional 
density ipip* in the space TL. Note that the equations (EH describe a classical 
field-Hamiltonian system on the space TL with respect to the symplectic form 
21m {p\ip), which determines Poisson brackets on the submanifold M C C: 

(2.8) {7 ,u} d (ipip*) = 2lm(H(ipip*)ip\A(ipip*)ip), 

where A(ipip*) = 8a{ipip*). 


2.3. A solution of equation (EH) can be written down implicitly for given 
initial condition ippto) = p by means of a time ordered mapping f5]: 


(2.9) 


ip(t) = exp <J —i PT(s,?/ ; (s)V’('S)*)ds| p, 


which generalizes the ordinary Feynman time ordered product. Thus, the 
solution of EH reduces in the given case to “disentangling” the expression 
(12.91) . i.e., to finding an operator VptHeuPP*) by means of which (12.91) can 
be represented explicitly as a function of p: 

( 2 . 10 ) ip(t) = V(t,t 0 ,pp*)p 


Sometimes, this can be readily done only for certain conditions p. For 
example, if for every v = p*p there exists a vector p = p v such that 


( 2 . 11 ) 


H(8,p„pt)p v = u(s,v)p v 


directly for all t, then the solution (l2~T0l) has the form 


( 2 . 12 ) 


V’(t) = exp|— i J co(s,v )d 


<Pv 


This fact is a special case of the following assertion. 
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Proposition 2. Suppose Equation. has one or several integrals 7 Tj(ifif*) 

of the form nj(ifif*) = if*Pjif, where Pj are self-adjoint commuting oper¬ 
ators in hi, and suppose there exists an extremal if = <p p , p = {pj} of the 
functional 7 (t,ifif*) which is the same for all t and satisfies the conditions 
TpPjTp = Pj> 

(2.13) H(t, p p p* p )p p = Y v’Pj'Ppi 

0 

where v 0 are parameters conjugate to pj (Lagrangian multipliers) determined 
by the conditions 

(2.14) u J (t., p) = dk g p P \ h(t,p) ='y(t,tp p tPp)- 

Then the solution of Equation. with the initial condition if(0) = p p 

has the form 

(2.15) if(t) = expj-i^ v(s,p)Pds| 
where vP = 'Yf Pj ■ 

This can be proved by direct substitution of (12^51) into m with al¬ 
lowance for the invariance of 7 under the transformations if i—> e iqp ^ for all 
real parameters q = The parameters q\ zP, p 3 have the meaning of 

generalized coordinates, velocities and momenta respectively. 

Definition 1. The solution 1051) of Equation. CT is called a generalized 
soliton with momenta pj. 

To the soliton there obviously corresponds the case of existence of 

an extremal p v for the integral ttq (if if*) = if*if generated by the identity 
operator Po = P 

Suppose, for example, the space TL is the space of functions if(x) on 
3 x with scalar product (p\if) = f ip (x) if (x) dx, and the functional 7 
is invariant with respect to the group of displacements: 

7 {t, V’qV’q) = 7 (t, V’V’*) W>q(x) = iffx -q),q£ M fZ . 

In this case, the conditions of proposition © are satisfied by the ordinary 
momenta 

-Kjfifif*) = -i J if(x)djif(x.) dx,j = 1 , ...,d, 

and one can expect that ordinary soliton solutions would exist - travelling 
waves with momenta p = {p 3 }: 
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This last holds, for example, for the ordinary (scalar) Hartree-Fock equation, 
for which the Hamilton operator has the form 


(2.16) 


= — d 2 + J cu(x — x') |t/i(x , )|" dx 7 , 


where d 2 is the Laplacian with respect to x £ M d , and where u(x — x') is the 
potential of the two-body interaction, for example, the Green’s function for 
the equation d 2 u(x) = 47 r ^(x)) 2 . Our treatment also encompasses more 
general cases of Equation. (EH), for example, if d 2 is a square with respect 
to some indefinite metric in An interesting case arises when M. d is four- 
dimensional space-time, for which d 2 is the d’Alembert operator □. 


2.4. In the general case, if there is a discrete basis in TL , it is convenient 
to identify the vectors ip* and tp with covariant and contravariant tensors 
of the first rank ip x = Tp fixing this basis. In this notation, the form 
i/)*tp = (ip |?/>) is simply the tensor contraction ifj x ip x such that the mapping 
-0 e-> ijj* is the operation of complex conjugation and lowering the super¬ 
script by means of an antilinear metric tensor J : ij: x = J xy ip y , = J xy i^ y . 

At the same time, the operators p E C, A e A are represented by ma¬ 
trices p x . A y , the bilinear form ( p,A ) is the contraction p y x A x , and the 

derivative H{p) = S'y(p) is given by the partial derivatives H x (p) = d ^ p J . 

Further, we shall encounter tensor powers p® n , 5® n , which in our notation 
are the products p y Xl ...p y Xri , 5 y ^...8 y ™ = *i dp —• Formally, the tensor 

notation can also be used even in the continuous index case if the partial 
derivatives are understood as variational derivatives and contraction as an 
integral with respect to a given measure; however the matrix elements are 
then represented by generalized functions. 


3. Uniformization and Second Quantization 

3.1. We shall say that a Hamiltonian system (EH is uniformizable if the 
Hamilton functional y(f, p) is a polynomial: 

N 

(3.i) ^t, P ) = Y,^(p® n , w(n) w) 

71=1 

or entire function (N = oo). Here, p® n are tensor powers of the element p 
G L generating the symmetrical tensor-spaces £„ 3 , W^ n \t) are elements of 
the dual algebras AA"> = C* n , and (p® n , are linear forms on AA') 3 

AA) equal to the products (p, A \)... (p, A n ) on the A ( - n )-generating elements 
AA) = Ai®...®= Af n . Note that WA) are represented as the symmetric 
tensors in any discrete basis, and are uniquely determined from 7 (p) as n-th 

^That is C n is the minimal Freshet subspace of the projective tensor power £® n con¬ 
taining all p®", p £ £. 
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derivatives S^j(p) at the point p = 0 defined as the tensor power 5^ = 5® n 
of the derivative 5 in m such that 

— e {p,A) A®n 

for any A € A. The symmetrical spaces *4^ like the full tensor power 
spaces A <z>n are Hamiltonian (Lie-Jordan) algebras with respect to products 
//(”) • A^ n \ { H ^ , A^ } defined on the primitive elements H <n ' 1 = (££)™ =1 Hi , 
A^ = <g)? =1 Ai as 

n n n 

(3.2) tf("> -A^ =(g)H i -A i , {h-("),hW} =Y / <g) H i-A i ®{H j ,A j } 

i= 1 7=1 2=1 

such that Hf n ■ Af n = {H^)® 11 , {H® n ,Af n j = nH® {n ~ l) • Af (n_1) 
{H^A,}. 

The symmetry of W^ n \t) makes it possible to write the Vlasov Hamil¬ 
tonian H{p) = S r j(p) determining the Vlasov equation (12.311 . in the form 


N—l 

(3.3) H {t,p)=Y J - ] {p® n M n+1 \t))i 

71=1 

where (p® n , the elements of the algebra A, are equal to (p,A\) n A 

for A^ n+l ^ = Ai (g) ... (g) A\. An example of such an expression with N = 2 
gives the operator TO if is replaced by an arbitrary density operator 

p = E-MhV’*- 

Definition 2. Let £ > 0 be a parameter; an e-uniformization of the classi¬ 
cal Hamiltonian system iEU with the Hamiltonian functional 7 ( t, p) is the 
system of the Hamiltonians 

(3.4) 

HW(t) = ^(1 + 5e)® n 'y(t, p) | p =o= (jji) £m ~ 1 ( l(n ~ m) 

m= 1 ' ' 

as the real elements in the Hamiltonian algebras A^. Here /( n_m ) g> S W^ 
is the symmetrization of the tensor products /( n_m ) <g> w( m ) 0 f the identities 
j(n—m) = / (g)... (g) / of the algebras A^ n ~ rn ^ and the elements 6 A^. 

For example, the uniformization of the Hamilton operators of the 

Hartree equation m leads to the system of operators 

n n i\ — 1 

(3.5) H (n ^ = + IZ w (xu -Xi 2 ), 

2 = 1 21 = 1 22 = 1 

where df are the Laplace operators with respect to the variables x* e R^. 
For four-dinrensional space-time R 4 , for which d 2 is the d’Alembert operator, 
the uniformized Hamiltonians are ultrahyperbolic operators for n ^ 2. 
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3.2. Note that the Hamiltonians cn determine in the spaces C n the linear 
equations 

(3.6) 

which are abstract Liouville (or von Neumann) equations understood as 
predual to the abstract Heisenberg equations 

( 3 . 7 ) = ^ H (n) {t ^ A n {t) y 

so that the relation l^p n (t), Ay 1 ^ = A^) (t)) holds for all initial = 

A (n) (t 0 ) G A (n) and p° n = p n (t 0 ) G C n . 

In order to develop a compact functional derivative representation of all 
these independent equations and find their relation to the abstract Vlasov 
equation we introduce the representing functionals 

(3.8) a(p) = exp {-) (p, /)} £ (p®”, A<">) , 

by means of which the elements AA) G AA> can be calculated as the deriva¬ 
tives 

(3.9) A("> = £ "- 1 ^exp|i( i 0 ,/)|a(p) | p=0 = ^(1 + e5)® n a(p) | p=0 . 
Here, A^ is a complex number A^ 0 ) G C = A® . 

Proposition 3. The system of equations \3. 7| ) can be expressed in the form 
of the pseudodifferential equation 

(3.10) 9aP - = ^ j p,H (t,p- ^1+ £5^ | ,5^ ot(t,p ), 
defined in terms of the symbol-Hamiltonian 

H (t, p ■ (I + eA)) = exp je (f>, 5 ■ A^ j H(t, p), 

where the superscripts indicate the order of application of the operators of 
(tensor) multiplication by p and the derivation with respect to p. The equa¬ 
tion (SUW has the limit at £ —> 0 of first order differential equation 

d 

Qjfxfo P) = ({P’ H (*> P)} > 5 ) «(*> P ) 

giving the weak form of the Vlasov equation for a (t, p) = (p ( t ), A). 

The proof of this proposition is contained in the derivation of the formula 
(expansions) (13.131) and (13.1 41) in the Appendix. Here we only note that the 
symbol 

G(t, p, A ) = <{p, H (t, p ■ (I + eA))} , A) 
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of the equation dm which is determined on the exponential functionals 
a (p) = exp (p, A) by 

(3.11) G (t, p, e {p ' A) = e {p ’ A) G (t, p, A ), 


depends analytically on e as a power series expansion in e. Comparing 
EH with El , we readily see that for e = 0 the derivative EH simply 
coincides with the derivative of a(t, p) = a(p(t)) along the trajectories of 
the Vlasov equation El- Expanding the operator G (t, p , 5) in powers 
of e and substituting in Equation. EH the solution in the form of the 
series a(t,p) = Yl'kLo £ka k(ti p), we obtain by the method of successive 
approximation the following system of coupled equations: 

(3.12) 


9ak{t) _ y' 1 

dt (n + 1)! 

n =0 




Qk—nify i 


where H( m \t,p) = Thus, the Vlasov equation (12.31) is a Hamil¬ 

tonian system of the characteristics for the zeroth-order approximation cto (t) 
with initial condition a (to) = «o- 


3.3. In general, the uniformization procedure can be regarded as a linear 
functor from the classical algebra of entire Hamilton functions 7 (p) with 
ordinary product 7 (p)a(p) and Poisson bracket (12.51) to a new Lie-Jordan 
(Hamiltonian) algebra 21 of these functions with not necessarily associative 
product 

(3.13) (7 -a)(p) = ^7 (p- + V(p) 

fi — ^ 

71=0 

and uniformized Poisson brackets 

(3.14) {7 ,«}(p) = (p,{MP' (/ + £$)) ,<Mp)}) 

= 7 (p),^ n «(p)}>- 

n=0 n ' 

The latter in accordance with El determines for a fixed Hamilton func¬ 
tional 7 (t, p) the uniformized Hamiltonian derivation Ell on the algebra 
21, in which H(t,p ) = <Py(t, p), since Equations (13.81) and (13.91) establish an 
isomorphism between the algebra 21 and the direct sum algebra 
of symmetrized tensor powers AA'i of the algebra of elementary observables 
A since the products EH and EH are simply representing functionals 
for the products HA) ■ aA ) ; {HA) : aA)} of components of the elements 
A = ®^ 0 AA ) and H = ®)^ 0 HA) (this last follows from the proof of 
Proposition CO)- In the limit ei->0, the Jordan product EH multiplied 
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by e goes over into the ordinary pointwise product 7 (p)a(p), and the Lie 
product (ETH in this limit obviously coincides with the classical Poisson 
brackets (E3). 

In the case when the Jordan product in A is associative, the uniformized 
product 4H35D also remains associative. This is the case when A is the clas¬ 
sical algebra of single particle observables A(q, p). Then the uniformization 
is equivalent to the transition from the classical field Hamiltonian system 
described by the Vlasov equation itself to the classical Hamiltonians 
where £ is the coupling constant. In this sense, the uniformization procedure 
is the inverse of the procedure of derivation of the Vlasov equation from the 
Bogolyubov equation [fr. One could naturally call this replacement of the 
classical field system by a system with variable particle number as quantiza¬ 
tion, however, we shall avoid this terminology here since by quantization one 
always understands transitions from a classical to a definitely nonclassical 
system, whereas the example considered here shows that the uniformized 
system may remain classical. 

The same interpretation remains valid in the quantum case when A is 
an algebra of operators in Hilbert single-particle space Ti. The uniformized 
system of a variable number of identical quantized particles is, in contrast 
to the original field classical system, no longer classical, and in this sense 
the described procedure can be regarded as a certain quantization, and it 
can be given a Heisenberg from with generating commutation relations 

(3.15) [(p,H),(p,A)] = (p,[H,A]), 

if one introduces an associative and noncommutative product (7 • a) ( p ) as 
the representing functional for the products HA) . aA) of operators in the 
tensor powers TiA) of the Hilbert space Ti. At the same time, the com¬ 
plex extension of A is a linear noncommutative associative algebra isomor¬ 
phic to the algebra of operators A e S) = ®))V 0 of decomposed form 
A = ! 3 AA ), where AA ) are symmetric tensor operators in TiA) that leave 
invariant the subspaces Ti±^ C TiA ) of the completely symmetric (TL 

as well as the completely antisymmetric (ri^^j tensors tjjA) g TiA). This 
last means that the elements of the algebra A can be interpreted as the 
observables of a second quantized system of identical bosons or fermions 
depending on whether or not they are regarded as operators A in Fock 
spaces Sj± = that are symmetric or antisymmetric under per¬ 

mutations of particles. Such Fock representations of the algebra A are ob¬ 
tained by means of the functionals (E3) in accordance with the formula 

A = )a , where a±, a)j_ are vector operators of annihilation and 

creation of bosons (with plus sign) and fermions (with minus sign) acting on 
the Fock spaces S) + , respectively, and the numbers 1 and 2 above the op¬ 
erators indicate the normal order, in accordance with which the monomials 




NONLINEAR HAMILTONIAN KINETIC EQUATIONS 


13 


^ ^aa*^ , are written as the normal products a*® n a ® n . The 

linear observables a ( p ) = (p , A) as functions on A 9 A are then represented 
in by operator-valued forms h(A) = a±Aa±, and the commutation rela¬ 
tions (EZE3) can be deduced directly from the canonical commutation rela¬ 
tions for the operators a±, a^.: 

(3.16) [h(H),n(A)\ = [a±Ha±,a±Aa±] = a?j_ [H, A] a± = h ([ H , A]). 


Hamiltonians H = ~ 7 ^eaa*^ of the nonlinear form EH) can be represented 
in the Fock spaces S)± by expansions in the coupling constant e: 

N n _i 

(3.17) H(t) = -— —a’^ >n W^ n ' > (t)a® n . 

n=l U - 


Thus, uniformization of the quantum Vlasov equation can be regarded as 
the algebraic basis for a simplified method of second quantization which 
avoids explicit introduction of unphysical - odd - observables (which do 
not commute with the operator of the total particle number n ) and 
takes into account in a unified manner for Bose and Fermi systems the 
principle of indistinguishability of identical particles. The calculus (13131 
of the representing functionals ESI together with Equation. EH 
combines the holomorphic calculus of the operators A = 0 A in the 
fock spaces S)± developed by Bargmann m for the case of Bose-Einstein 
statistics and Berezin EH for Fermi-Dirac statistics; at the same time, the 
commutation relations ( 13.1511 replace the canonical commutation relations 
for the operators a ±, alj_. 

All that we have said remains true in the more general case when the 
metric of the space Tt can be indefinite except for the interpretation of uni¬ 
formization as second quantization since the spaces Sj = like their 

subspaces S)± need not be Hilbert spaces, but only, like Tt , pseudo-Hilbert 
spaces. The representations ® of the algebra A by the operators ( 13.81 
in the spaces TL^) are then representations with indefinite metric, and are 
not Fock representations. For this reason, the dual elements p = <S) n p n , 
where p n E C n are operators of the form p n = )>V cannot, in 

general, be interpreted as statistical operators of n-particle states even if 
Aj ^ 0 and Aj = 1. Nevertheless, investigation of the linear equations 
EH or the (abstract) Schrodinger equations 


(3.18) 


1 dt 




1 dt 




into which they decompose for initial conditions of the one-dimensional form 
p n (t 0 ) = enables one to extract the entire information about the 

dynamics of the uniformized system. It is sufficient to consider only the 
symmetric or antisymmetric n-tensors E Tt^ n \ for which the products 
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v> (n V (n) * are elements of the space C n of symmetric n-tensors p n since 
equations (ETCH) preserve the symmetry of ip^ by virtue of the fact that 
the Hamiltonians H^ G commute with the permutation operators. 

4. Uniformized Solutions of the Hartree Equation 

4.1. Here, we consider a method based on uniformization for disentan¬ 
gling the Time ordered mapping uni, this giving solutions of the Hartree 
equation EH) with Hamiltonian of the form EH, where p = ipip* and 
{p® n , = ip*® n w( n+1 ' > ip® n . To be specific, we shall assume through¬ 

out that the forms ip*® n j\f n ) are c-numbers, although they could also 
be regarded as even elements of a Grassmann algebra with ip and ip* re¬ 
garded as anticommuting generators. At the same time, ip*® n ip® 71 = 
ip*A\ip ■ ■ ■ ip*A n ip for AAi, and ip*Aip = (Ap\ip). 

We shall regard the differences H^ n+1 ' l (t) — H^ n \t) <g) I as operators in 
TL^ <g> TL and denoted by £}( n+1 \t, to) the operators in TL^ <g> Tt obtained 
by there transformation: 

(4.1) rt n+1 \t,to) = £/ (n) (Mo)* -H<- n \t) ®l) U in \t,to), 


where U( n \t,to) are unitary operators in n" ' 1 determining the solutions of 
equations EH: 


(4.2) 


U^\t,to) = 


exp < —i 


r* s (n) 

H (s)ds 


>t 0 


Definition 3. We give the name e solution of equation with the op¬ 

erator eh and initial condition ip (to) = ip to the following expression 
ip(t ) = V(t, to, ipip*)ip, where 

f 1 ) 00 p n—l 

(4.3) V(t,to,<P<P*)<P = exp l --ip*ip \ ^ —^~<p*® n V ( - n+1 ' ) (t,to)v® n , 

' ' n =0 

and V^ n+1 \t,t o) are unitary operators in TL^ ® hi with generators 
i.e., the time-ordered products 


(4.4) 


V^ n+1 \t,t 0 ) = 


exp < —i 


'to 


(n+l) 

H (s,to)ds 


The justification for this definition is the following proposition, which 
reduces the problem of disentangling the Time-ordered mapping (12.91) for 
uniformized systems to calculation of Time-ordered products EH- 

Proposition 4. The solution of the uniformized equation \2. 7| ) with any ini¬ 
tial condition ip(to) = p coincides with the £ —> 0 limit of the corresponding 
e solution. 

This proposition is readily proved by noting that the Time-ordered prod¬ 
ucts EH are found by disentangling the operator (jjA) <§> I) with the 
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generator (g) I from the operator [/( n + 1 ) with generator H^ n+l ^ = 

<g)I+^( n+1 ), 

(4.5) V^ n+1 \t,t 0 )= (t/W(Mo)®l)”V (n+1) (Mo). 

This, the e solution essentially reduces to finding the unitary operators 
representing the solutions of the Schrodinger equation in the form 

(t,to)ip( n ^ for any initial condition p € TL^. The operators 
m multiplied by y/n are none other than the nonzero matrix elements of 
the vector operators of annihilation of bosons a(t ) in the Heisenberg rep¬ 
resentation, and the e solution as a function of p, p* is the symbol of the 
vector operator sjea(t) = ip(t) ordered normally with respect to the initial 
operators p = y/ea + (t 0 ). For e = 0, the initial vector operators p coincide 
with the numerical vectors p, and the e solutions coincide with the solutions 
of the Hartree equation since the Heisenberg equation for the symbol ip(t) 
in the representation of the instantaneous time t coincides exactly with the 
Hartree equation. 

4.2. The disentangling of the Time-ordered mappings in this way requires 
the solution of all the Schrodinger equations (EHsl) for n = 1,2, .... It is of 
interest to ask whether one can find not all but at least some solutions of 
the Hartree equation by solving not all but, for example, two neighboring 
equations (13T81) . 

Let p^ n \ p^ n+l ^ be (generalized) stationary solutions of Equations. (15.1811 
corresponding to the eigenvalues X n (t), A n+ i(t) of the uniformized Hamil¬ 
tonians H^ n \t), H^ n+1 \t). One can evidently write the matrix element 
ip(t) = <^( n )*lA n + 1 )(t, t 0 )<^ n+1 ) of the operator (14.51) in the form (EH, 
where 

(4.6) p( n >p( n+1 '> = v = en, uj(t, v) = X n+1 (t) - X n (t). 

If all quantities have a limit as e —> 0 , then one can expect that in 
this limit the given matrix element defines a certain stationary solution of 
the Hartree equation, while for e ^ 0 it gives an approximate solution of 
this equation. 

Generally speaking, suppose the Hartree equation JO) has one or several 
integrals of quadratic form p*Pjp, including Po = I, and suppose |p) are 
(generalized) eigenvectors of the operators ea*Pja and the Hamilton opera- 

/ i 2 *\ / \ 

tor 7 I saa ) in the space f) = 0 corresponding to the eigenvalues pj 

and h(t, p). We denote the integrals of the matrix elements of the Heisenberg 
vector operators a(t ) by 


(4.7) 


(p|a(t) |p ; ) dp'. 
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Proposition 5. The vectors ra can be expressed in terms of the initial 
vectors ip p = if (to) in accordance with 

(4.8) ip p (t) = expj-^ (h{s,p + eP) - /i(s,p))ds| ip p . 


Proof. For any Heisenberg operator Aft) we introduce in the space Sj the 
function 


(4.9) a ft, p, q) = J exp | ^-q- (p - p') | (p| Aft) |p') dp'. 

It is easy to see that the Heisenberg equations ezd lead to the following 
equation for the function (El: 


(4.10) 


.da(t,p,q) 1 


dt 


d 


= -(h(t,p + ie— ) + h (t, p) ) a ft, p, q) 


In particular, for the vector function 


(4.11) V’p.q(t) = J exp | ^q- (p - p') | (p| aft) |p') dp' 

Equation. EH can be written in the form 
dib n Jt) 1 

(4.12) i -- = - (h ft, p + eP) - h ft, p)) ^ p>q (t), 

where we have noted that i = Pjif pq ft) in accordance with the re¬ 
lation [a*Pja. a] = Pjd. The solution of the linear equation EH with the 
initial condition tp = V’ P ,o(^o) can be written in the form (14.81) . as we wish 
to prove. | 


If the function h ( t , p) in the limit e —> 0 has a p-differentiable limit, then 
El obviously takes the form of the soliton solution EH of the Hartree 
equation. In this case, for e A 0, the vector function El can be called an 
e soliton, or an extremon ,3). 


5. Appendix 

Taking into account equations El and El. we write the representing 
functions 


(5.1) (7 -a)(p) 

(5.2) {7,«}(p) 


k J n =0 

exp {-I <„./>} E^( P ®", {«<’),A«}) 

v ' n=0 


for the products ■ A^ n \ {H^ n \ A^} in the symbolic form 




E 

n=0 


1 


n\e n+1 


(p^Al + eSA 6 


(/ + e6 2 ) (Z,n ) 7 {pA a (p 2 ) | Pl=0 = P2 , 
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oo _. 

+ + Ipj-o-pi , 

n=0 

where (5i and (52 are the derivatives with respect to p 1 and p 2 respectively. 
Using further Equations. Id. 21 for the powers = (I + £(5i) 0n , A ^ = 

(/ + we obtain 

(5.3) (p® n ,(I + e8 1 f n -(I + e5 2 f n ) = 

(P, (I + e5 1 )-(I + e5 2 )) n , 

(5.4) (p m ,{(I + e<5i) 0 ” , (.I + £<5 2 ) 0 ”}} = 

n (p, (I + edi) • (/ + £( 5 2 ))” _1 (p, {e<5i, e<5 2 }). 

Substituting this into the preceding formulas, we obtain 

(7 • a) (p) = i e l<oV+rf.)-a+«&)-') 7 (ft) „ (p 2 ) | pi=0 , ft , 

{ 7 , a} (p) = e 7 (pXi+^iHi+eS 2 )-i) ^ | Ji7 ( Pi ) 5 s 2 a (p 2 )}) |p 1 =o=p 2 . 

Multiplying the exponentials and noting that ^/3(p') | p '=o= P(p), w e 

obtain 

(5-5) (7 • «) (p) = ^e £{p ’ Sl 52) -f {p^ a (p 2 ) \ Pl = P =p 2 

(5.6) { 7 , a} (p) = e £{p ’ Sl S2) (p, {< 5 i 7 (p x ), 5 2 a (p 2 )}) I Pi =p=p 2 ■ 

Expanding the exponential Q in a series in e, we obtain the expansion 
dnnsi)- The first formula (EH, from which Proposition IB follows, is ob¬ 
tained by applying the second formula (15.51) to (15.61) . and the expansion 
El in powers of e can be readily obtained by using formulas of the type 
(tPli and (fOI) in reverse order in the exponential (EH). 
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